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ABSTRACT 

We have studied superionization and X-ray line formation in the spectra of £ Pup using our new 
stellar atmosphere code (XCMFGEN) that can be used to simultaneously analyze optical, UV, and 
X-ray observations. Here, we present results on the formation of the O VI AA1032, 1038 doublet. 
Our simulations, supported by simple theoretical calculations, show that clumped wind models that 
assume void in the interclump space cannot reproduce the observed O VI profiles. However, enough 
O VI can be produced if the voids are filled by a low density gas. The recombination of O VI is very 
efficient in the dense material but in the tenuous interclump region an observable amount of O VI 
can be maintained. We also find that different UV resonance lines are sensitive to different density 
regimes in £ Pup : C IV is almost exclusively formed within the densest regions, while the majority of 
O VI resides between clumps. N V is an intermediate case, with contributions from both the tenuous 
gas and clumps. 

Subject headings: stars:early-type - stars:winds, outflow - X-rays: stars - X-rays:individual:£ Pup 



1. INTRODUCTION 

One of the surprising discoveries of the Copernicus 
satellite was the strong P-Cygnii profiles of superions, 
such as O VI and N V, in the FUV spectra of many O 
and B stars (Snow & Morton 1976). The only viable ex- 
planation for the presence of O VI is Auger ionization 
by X-rays from O IV (Cassinelli & Olson 1979) which 
is the dominant form of oxygen in many O-type stars. 
The X-ray emission, necessary for Auger ionization, was 
later detected by the first X-ray telescopes (e.g., Harnden 
et al. 1979; Seward et al. 1979). 

The origin of the stellar X-ray emission was another 
enigma until the "wind-shock" mechanism (Lucy & 
White 1980) became the accepted explanation. Massive 
stars posses strong line-driven winds in which the ma- 
terial is accelerated by numerous C, N, O, and Fe tran- 
sitions (see e.g., Pauldrach & Puis 1990; Castor et al. 
1975). It was known from the conception of the line- 
driven wind theory that such flows are unstable and 
prone to the formation of dense clumps and shocks (see 
e.g., Owocki & Rybicki 1991; Lucy & White 1980). The 
large-scale flow energy is converted to heat in the shock 
fronts producing high temperature plasma. Numerical 
simulations confirm (Feldmeier et al. 1997; Owocki et al. 
1988) that at least the soft X-ray emission of early-type 
stars can be explained by this mechanism. 

Evidence for density inhomogeneities (or clumped 
winds) is provided by variability studies of both WR 
stars (Lepine & Moffat 1999, and references therein), 
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and O stars (e.g., Eversberg et al. 1998; Lepine & Mof- 
fat 2008). Further, density inhomogeneities allow the 
electron-scattering wings of emission lines to be reduced 
to the observed level while maintaining the strength of 
emission lines (Hillier 1991; Hamann & Koesterke 1998; 
Hillier & Miller 1999). More recently, Crowther et al. 
(2002) and Hillier et al. (2003) found that they could 
not simultaneously fit the Ha and P V AA1120 profiles 
in normal O supergiants without assuming an inhomo- 
geneous density distribution in the wind. Using a more 
statistical approach, Massa et al. (2003) showed that the 
phosphorus ionization structure was consistent with ex- 
pectations only if lower than conventional mass-loss rates 
were used in their analysis of P V AA1120. Additional ob- 
servational evidence for wind clumping come from Chan- 
dra and XMM-Newton high-resolution X-ray spectra of 
O stars. These spectra revealed that X-ray lines suffer 
less absorption in the wind than predicted by "smooth" 
models (e.g., Kramer et al. 2003; Cohen et al. 2006; Wal- 
dron & Cassinelli 2007, and references therein). 

Superionization has received only limited attention 
since the work of MacFarlane et al. (1993; 1994). The ef- 
fect was introduced into modern stellar atmosphere codes 
(e.g., WM-basic, Pauldrach et al. 2001), but we are un- 
aware of any work that has revisited the question in light 
of the high-resolution X-ray observations, improved X- 
ray emission calculations and the new results on clump- 
ing. With improvements to CMFGEN, we are developing 
tools and techniques to move towards this goal. As part 
of this effort we discovered that the interclump medium 
is crucial to explain the O VI doublet profile in £ Pup. 
In this paper we demonstrate the effect and discuss its 
implications. In §2 we briefly describe our code, the ob- 
servations we used, and our models. We present and 
discuss our results in §3— §4. 

2. OBSERVATIONS, TOOLS, AND MODELS 

We coadded Copernicus Ul scans for £ Pup to create 
the observed profile in Figure 1. All but two of these 
scans were blocked Ul scans that minimize stray-light 
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as described in Rogerson et al. (1973a, b). Our data re- 
duction (see, Zsargo & Hillier 2008) produced an O VI 
profile that strongly resemble those in MacFarlane et al. 
(1993) and Morton (1976). 

We have used XCMFGEN (Zsargo & Hillier 2008) to 
solve for the ionization balance and the non-LTE level 
populations in our stellar models. XCMFGEN is a new 
version of CMFGEN (Hillier & Miller 1998) that can 
perform X-ray emission calculations in addition to the 
original CMFGEN tasks. We use APEC (Smith et al. 
2001) and its accompanying database of X-ray cross sec- 
tions (APED) to calculate non-LTE level populations 
and emissivities in the X-ray emitting plasma. 

Bouret et al. (2005, 2008) used CMFGEN to derive 
stellar parameters for seven O supergiants. We used their 
values for £ Pup to construct the models listed in Table 1. 
Our improved models specifically addressed formation of 
O VI and X-ray lines (Zsargo & Hillier 2008). The X-ray 
emitting plasma was distributed in the wind and its cmis- 
sivity was characterized by three plasma temperatures. 
The parametrization included a volume filling factor (fx) 
and an initial turn-on radius (Ro), as in Owocki & Cohen 
(2001). The volume filling factor was constant beyond 
Ro and was adjusted until the prescribed Lx/Lbol was 
met (see Table 1). Note that Lx is the X-ray luminosity 
attenuated by the wind and not the intrinsic luminosity. 

To treat clumping we follow the "volume filling factor" 
approach where the material is compressed into a fraction 
of the available volume with void in between (see e.g., 
Hillier & Miller 1999; Hamann & Koesterke 1998). The 
ratio of the volume filled with material (clumps) to the 
total is 



fdir) = foo + (1 - foe) ■ exp(-v(r)/vd) 



(1) 



where foe and v c i are free parameters. Eq. 1 is an ad-hoc 
formula, motivated by hydrodynamical simulations, and 
provides a smooth wind at low v(r) velocities (r ~ i?*). 

In our models, fx and f c i are two independent pa- 
rameters despite the fact that both are referred to as 
"filling factors" . The two should be related since both 
of them are the results of the wind instability. However, 
our understanding of the line-driven winds is too poor to 
formulate this relationship. 

3. RESULTS AND DISCUSSION 

Figure 1 shows our model calculations together with 
the Copernicus spectrum. Auger ionization produces 
strong O VI features in the case of the smooth wind 
model, albeit a bit weaker than observed. A better fit can 
be achieved by increasing the oxygen abundance and by 
assuming larger turbulent velocities near . A detailed 
analysis of £ Pup and the simultaneous fit to the X-ray 
and UV spectra, will be presented in Zsargo & Hillier 
(2008). The most important feature of Fig. 1 is the ab- 
sence of an O VI profile in the clumped wind model. This 
cannot be rectified by fiddling with, for example, the oxy- 
gen abundance. Adjustments at the order of magnitude 
level would be necessary to produce a visible O VI pro- 
file; this indicates that something fundamental is wrong 
with the simplified "clumped" model. 

There is a simple theoretical explanation for the weak 
O VI P-Cygnii profiles in clumped models. Since the 
O VI AA1032, 1038 doublet forms only in the wind by 
redistribution of the stellar radiation, the optical depth 



TABLE 1 

Model Parameters (representative for C Pup). 



Photosphere /Wind 

R» a 



log g (cgs) 

v turb 

v sin i 

M 

Joo 
V C ! 

Plasma 

Ro c 

Lx/Lsol 



L pi 
[pi 



i a 



19 R 
39,000 K 
3.6 

7.0xl0 5 L 
12 km s" 1 
240 km s" 1 
2300 km s" 1 
0.9 

Smooth Wind Clumped Wind 

7.6X10" 6 Moyr" 1 1.7xl0" 6 M yr" 
0.05 
150 km s" 1 



1.5 R» 
1.5xl0~ 7 
4.6xl0 6 K 
2.4xl0 6 K 
1.7x10 s K 



Note. — Stellar parameters arc from Bouret et al. (2008), unless 
noted otherwise. Plasma parameters are from Zsargo & Hillier 
(2008), unless noted otherwise. For simplicity, the same emission 
measure EMx = J n e njjdV was assumed for each T p ;. 

a Value is from Hillier et al. (1993). 

b Power of the CAK velocity law (Castor et al. 1975). 

c Value is from Kramer et al. (2003); Leutenegger et al. (2006). 




Fig. 1. — Comparison of the Copernicus observation of O VI 
AA1035 (rugged line) with those calculated by XCMFGEN for a 
smooth wind (red), for a clumped wind model (blue), and for the 
interclump medium (green). See Table 1 and §3 for details. 

controls the shape and strength of the profiles (see e.g., 
Lamcrs et al. 1987). In a spherical and accelerating flow 
it is appropriate to write the optical depth as 
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where the atomic parameters have the usual meanings 
and arccos(/Lt) is the angle between a line-of-sight and 
the radial direction. The average O VI density can be 
written as 



no vi = A n q Q vi ~ A ( 
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by using the oxygen abundance Ao, the O VI fraction 
Qo vi, and the average particle mass m p . Note, that 



On the importance of 



no vi is the mean ion density over a Sobolev length, so 
no vi = no vi, d fd = A n c i f c i q vi (4) 

must be used if the wind material occupies only a fraction 
fd of the available volume with n c i number density. 

The key unknown in the expression for to vi is the 
O VI fraction. We can estimate qo vi by a simplified 
rate equation: 

dn Q vi 



dt 



no iv a(Jx) - n vi n e 7 = 0, 



(5) 



where a(Jx) is the X-ray-flux-dependent rate for Auger 
ionization, and 7 is the effective recombination rate co- 
efficient. Solving for the ion fraction, 



no vi _ qo vi _ ct(Jx) 



(6) 



no iv qo iv n e j 
Under the approximation that qo iv ~ 1 and n e ~ n, we 
find that <? vi = e ^~p-, or 



f»o vi = A n q vi 
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(7) 



The ion density of O VI, therefore, is independent of the 
wind density, and thus the mass-loss rate. Substituting 
this in the expression for the optical depth, we find for a 
smooth wind 



TQ VI, sm 
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r \ dr 
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The Sobolev optical depth in a smooth wind is thus also 
independent of the mass-loss rate. 

Let us generalize this to a clumped wind. Equations 6 
and 7 still give the O VI fraction and density, but these 
are the values for the dense clumps, and not means over 
a Sobolev length. Thus we need to use Eqs. 2, 4, and 6, 
together with the assumptions of qo iv ~ 1 an d n e ~ n 
to derive 

TO VI, cZ ~ fcl t O VI, sm (9) 

in a clumped wind. 

It is clear now why clumped models fail to produce 
observable O VI lines. The volume filling factor for £ Pup 
{fd ~ 0.05 at r » J?*) decreases to vi.cZ well below 
unity. The results of our XCMFGEN simulations fully 
support the above calculation. Figure 2 shows the radial 
optical depth of the important UV resonance lines in our 
models (Table 1). The curves show what is predicted 
by Eq. 9 for small f c f, the O VI optical depth is nearly 
two orders of magnitude lower than in the smooth wind 
model. The dense clumps, therefore, contribute little or 
nothing to the O VI profiles. 

There are two other factors that further decrease the 
O VI optical depth, and hence adversely influence the 
model O VI profile. First, Bouret et al. (2008) found 
sub-solar oxygen abundances for all stars in their sample 
(a factor of five difference for £ Pup); this affects only 
comparisons with earlier calculations. Second, the X-ray 
flux available for Auger ionization tends to decrease when 
fd is lowered; the wind absorption is smaller, hence less 
X-ray emission is needed to meet the observed Lx /^bol- 

What are the implications of our failure to produce 
observable O VI lines with clumped winds? Is the wind 
smoother than suggested by Bouret et al. (2008)? Their 




Fig. 2. — The radial optical depth as a function of normalized 
wind velocity for the stronger component of the N V AA1240 (red), 
C IV AA1550 (blue), and O VI AA1035 (green) doublets. The 
solid and dashed lines are for smooth and clumped wind models 
(Table 1), respectively. 



filling factor for £ Pup is not extremely low; it is very 
similar to those found for other O supergiants (Hillicr 
et al. 2003; Massa et al. 2003). Further, increasing f c i to 
0.1 or 0.2 would not solve the problem. It is more likely 
that the use of the traditional "filling factor" approach 
in XCMFGEN causes the failure to produce observable 
O VI lines. Equations 8 and 9 suggest that the part of 
the wind with the largest filling factor will be the most 
important at producing O VI optical depth. Thus we 
need to consider the influence of the tenuous intcrclump 
medium (ICM) on the O VI profile. 

This is a non-trivial exercise; at least three compo- 
nents (dense clumps, ICM, and the hot plasma) need to 
be treated simultaneously for a fully self-consistent so- 
lution. Fortunately, the radiation field at any radius is 
almost independent of the ICM. Therefore, we can use 
J U (R) from our clumped model and solve the statistical 
equilibrium equations in the ICM only. We created such 
a model from the clumped model of Table 1. We scaled 
down the densities by a factor of /B to simulate the ICM 
at r >> i?*, and imported the radiation field from the 
original clumped model. We also assumed that the new 
model is smooth, reflecting the high volume filling frac- 
tion of the ICM. Effectively, we took the ICM component 
out of the real (clump+ICM+hot plasma) wind and built 
a stand-alone model for it. Since we imported the radia- 
tion field from the clumped model, the coupling between 
the ICM and the rest of the wind has been taken into 
account at least in the first order. 

The choice of density in the ICM is somewhat arbi- 
trary, but in the absence of reliable hydrodynamical pre- 
dictions it is still reasonable. There is a density contrast 
of four hundred between the dense clumps and the ICM 
at r > > i?* . This means that the ICM contributes little 
to the overall mass-loss rates. If we combined the mod- 
els for the ICM and clumps the total M would be nearly 
identical (within 5%) to that of the clumped model. 

The result of this experiment is displayed in Figs. 1 
and 3. It is obvious from both figures that the ICM has 
substantial O VI optical depth and contributes a strong 
P-Cygnii profile. Our results also offer a glimpse on the 
behavior of other lines. The C IV optical depth in the 
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Fig. 3. — Same as figure 2 but the solid curves are now for the 
interclump medium. See §3 for details. 

ICM is very low and thus most of the C IV AA1550 lines 
are produced by the clumps. The behavior of N V AA1240 
is very interesting. It appears from Fig. 3 that both the 
ICM and the dense medium contribute. 

Our results warrant an investigation of the physical 
conditions in the ICM and also its role in the formation 
of UV resonance lines in different environments (e.g., low 
and high density winds). For example, what is the den- 
sity contrast between the ICM and clumps? At very low 
densities the average charge state of oxygen (and other 
species) may be drastically different from the one in the 
dense regions. Also, how closely does the ICM follow the 
velocity of the clumps? Do O VI and C IV follow a dif- 
ferent velocity structure? Further, how are the ICM and 
hot plasma related? 

4. CONCLUSION 

In this letter we present our first results on superion- 
ization in clumped winds, and showcase the potential of 



the interclump medium to produce observable features. 
Clumped wind models that use the classical "volume fill- 
ing factor" approach (clumps with voids in between) can- 
not reproduce the observed O VI profile in £ Pup. The 
recombination of O VI is too efficient and the necessary 
fractional abundance cannot be sustained in the clumps. 
However, a tenuous interclump medium can contribute 
enough O VI to produce an observable O VI profile. Only 
a small amount of mass is necessary in the interclump 
medium, so its overall effect on the derived mass-loss 
rates is negligible. 

Our result highlights the need for improved treatment 
of clumping in the winds of massive stars. It is impos- 
sible to achieve a simultaneous fit to all UV P-Cygnii 
profiles with a single component wind model for £ Pup. 
Our simulations suggest that in £ Pup, different UV res- 
onance lines probe different density regimes. C IV is 
formed almost exclusively in the dense material, while 
O VI likely originates from the interclump medium. N V 
is an intermediate case with similar contributions from 
both components. In cooler O stars, when N 2+ become 
the dominant ionization stage, we might expect that N V 
shows the same behavior as O VI in the hotter stars. Ob- 
viously other possible effects of the interclump medium, 
in both O and W-R stars, should be investigated. 
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